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Abstract 

Chen's lemma on iterated integrals implies that certain identities 
involving multiple integrals, such as the de Bruijn and Wick formu- 
las, amount to combinatorial identities for Pfafhans and hafnians in 
shuffle algebras. We provide direct algebraic proofs of such shuffle 
identities, and obtain various generalizations. We also discuss some 
Pfaffian identities due to Sundquist and Ishikawa-Wakayama, and a 
Cauchy formula for anticommutative symmetric functions. Finally, 
we extend some of the previous considerations to hyperpfaffians and 
hyper hafnians. 

1 Introduction 

The free associative algebra over an alphabet A is naturally endowed 

with a commutative operation, the shuffle product LU which can be recursively 
denned by 

auilibv = a(uillbv) + b(aulUv) 

where a, b G A and u, v G A*. 

This operation has at least two interesting interpretations. First, it is 
dual to the comultiplication which admits the Lie polynomials as primitive 
elements (Friedrich's theorem, see, e.g., [Q). Second, it describes the mul- 
tiplication of iterated integrals (Chen's lemma 0). It is this last property 
which will be of interest to us. To state it in a convenient way, let us fix 
some vector space 7i of suitably integrable functions over an interval (a, b). 
For /i, . . . , f n G TC, we identify the tensor product F = fx ® • • • <g) f n with 
the function of n variables F(xi, . . . , x n ) = /^(xi)/^^) • • • fn(x n ). Let ( ) 
be the linear form defined on each 7i® n by 

(fl ® • • • ® fn) = / dx l / dx 2 '-- dx n fi(Xi)f 2 (x 2 ) ■ ■ ■ fnM ■ (1) 
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We now choose a family of functions (0i)igj, labelled by some alphabet J, 
and for a word u> = i\i 2 ■ ■ -i n G I*, we set 

(«;) = <g> i2 (8 • • • <g> 0i„) . (2) 

Then, Chen's lemma states that 

(u)(v) = (uUlv) . (3) 

Iterated integrals occur for instance in the power series solutions of lin- 
ear differential equations with variable coefficients. An example of current 
interest is provided by the multiple zeta values p0 |, and an interesting ap- 



plication of the shuffle product appears in ||, where a conjecture of Zagier 
|20| is proved by reduction to a simple shuffle identity. 

In certain cases, iterated integrals can be evaluated in closed form by 
means of de Bruijn's formulas M 



j ' ■■■ J det (fafe)) dxx--- dx 2n = Pf ( p ij)i< itj < 2n ( 4 ) 

where Pf (P) denotes the Pfaffian of the antisymmetric matrix P, 
Pij = J J [<Pi(x)4>j(y) - 4>j(x)(p i ( y y)]dxdy 

a<x<y<b 

and 



det (^i(xj)\ipi(xj)) dxx--- dx n = Pf (<5ij) 1 < iJ < 2n (5) 

a<xi<---<x n <b 

where Qij = J [<$>i(x)ij)j(x) — ipj(x)ipi(x)]dx and (<f>i(xj)\ipi(xj)) denotes the 
matrix whose zth row is [(f)i(xi), if>i(xi), <j>i(x 2 ), ipi(x 2 ), ■ ■ ■ ,<j>i( x n),ipi(x n )}. 
There is also a version of formula ([|) for determinants of odd order, 

J "J det (4>i(xj)) dx 1 ---dx 2n+1 

2n+l „i) 

V(_1)P+1 / <j) p (x)dxPf (P l3 )l<^<2n . 

~! Ja 



a<x\<---<X2n+i<b 

2n+l 



p=l 



The foundations of the theory of random matrices can be developed from 
these identities |Tj| . 
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2 The Wick formulas 

The above formulas are reminescent of identities which are usually called 
Wick formulas in the physics literature. These can be stated in various 
ways, and come in two flavors: bosonic and fermionic. 

The bosonic version amounts to a familar fact about moments of Gaussian 
measures. Let d/i(x) = C e~? xAx dx be a Gaussian probability mesure on WL N . 
We set (/) = f RN f(x)dfj,(x). Let /i, ...,/*.. . be linear forms. Then, the 
bosonic Wick formula reads 

it /• \ _ / if n is odd , . 

{h ■ ■ ■ tn) - | Hf {{lfi)) . .. h if n ^ even ? 1 > • 

the hafnian of a symmetric matrix of even order n = 2r being defined by 

Hf (A) = a u h a i2j2 ■ ■ ■ a lrjr (7) 
i,J 

where the sum runs over all decompositions of {1,2, . . . , 2r} into disjoints 
subsets / = {ik}, J = {jk}, such that i\ < • • • < i r , j± < • • • < j r and for each 
k, ik < jk- Replacing A by an antisymmetric matrix and multiplying each 
term of the sum by the signature of the permutation ixj\ ■ ■ -i r j r , one would 
obtain the Pfaflian of A. Actually, hafnians were introduced by Caianiello 
in order to emphasize the similarities between the bosonic and fermionic 
versions of the Wick formula. 

These similarities become particularly obvious when the fermionic Wick 
formula is stated in terms of Berezin integrals (integration over Grassmann 
variables 0). 

The Berezin integral is essentially a convenient notation for computing 
in Grassmann algebras. Let r\ = {t]i}i<i<N be a finite set of anticommuting 
symbols (i.e. r/ir/j +rjjrn = 0). We denote by /\ K (r)) the K-algebra generated 
by the rji. The Grassmann integral of an anticommutative polynomial / in 
the r\i is defined by 

<9 d 

dr]n--- drjifiVi, ■ ■ ■ , Vn) = -r i^—f(Vi, • • • , Vn) (8) 

where -£p acts on the Grassmann algebra as a left derivation (we push rji 
to the left, with a sign, and erase it). Therefore, in the fermionic calculus, 

integration and derivation are one and the same operation. 

it 

Let dfi(n) = Ce? vQv be a "fermionic Gaussian probability mesure" (Q 
being an antisymetric and non degenerate matrix, and j d/i(rj) = 1) and let 
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/i, be linear forms in the r/j. We define 

(/> = / dKv)f(v) (9) 
The fermionic Wick formula reads p| 

</l-/n> = ( L„ffV> -f niSOdd ( 10 ) 

W 7 1 Pf ((fifj})i<i,j<n ifniseven. v ; 

As we shall see, all these formulas are closely related, and reflect simple 
combinatorial identities in the shuffle algebra. 

3 Pfaffians and Hafnians in shuffle algebras 

3.1 Wick formulas in shuffle algebras 

The Wick formula (ffO) is equivalent to the equality 

E E ^•••/^■••^-n ( n ) 

n>0 ii<---<i n i \ i<j / 

More generally, to an antisymmetric matrix Q of order N over a com- 
mutative ring K, we can associate the following element of the Grassmann 
algebra f\ K {rf) over K 



T= n( i+ E^^) 

i \ i<j / 



(12) 

and this product being commutative we can express it as an exponential 

T = exp {I'vQv} = E "1 ( 13 ) 

where r\ denotes here the column vector (r]i). The coefficient (r]i\T) of r\i = 
Vh ' ' ' Vi2n i n ^ being ^ dr) i2n ■ ■ ■ dri^T, the fermionic Wick formula gives 

( VI \T) = Pf (Qi^k^n ( 14 ) 

Let us suppose now that K = R(A) Ui is the shuffle algebra over some 
commutative ring R. We will denote by Pf w(M) (resp. Hf w(M)) the Pfaf- 
fian (resp. hafnian) of a matrix M e R(A) Ui . 



Pfaffian and Hafnian Identities in ShufRe Algebras 



5 



Let (df) and (6$) be two sequences of letters of A, and set Q t j = dibj — ajbi. 
Remarking that 

(a h b h --- a ik b jk UJ a ik+1 b jk+1 ) n h n h ■ ■ ■ n ik+1 n jk+1 

ji,- Jk+i 

= XI a «i^ii ' ' ' a ii-i^ji~i a ik+i^jk+i a ifijl ' ' ' ^k^jkVhVji ' ' ' Vik+lVjk+i 

h,— ,3k+i l 

(15) 

since the left-hand side is equal to the right-hand side, plus a sum of terms 
which are symmetric in one pair of indices (j r ,j s ) an d are therefore zero, we 
find 

Yl e ( a ) a Mi) 6 M2) ' ' ' a i^k-x)K^)Vii ■ ■ ■ Vi2k 

= II ( 1 + ^3 ~ a ^ I ( 16 ) 

where e(er) is the signature of the permutation a, and formula flTl| ) leads to 
/] e(cr)a CT (i)& CT ( 2 ) ■ ■ ■ a a (2n-i)ba{2n) = Pf uj {dibj — %&i)x<jj<2„ (17) 

Example 3.1 With N — 4 we obtain 

01620364 — 01620463 — 01630264 + 01640462 + 01640263 - 01640362 
—02610364 + 02040263 + 02630164 — 02630461 H — 02640164 + 02640361 
—03620164 + 03620461 + 03610264 — 03610462 — 03640162 + 03640261 
—04620361 + 04620163 + 04630261 — 04630162 — 04610262 + 04610263 
= (^1^2 — ^261) LU (0363 — 0464) — (0163 — a 3 6i) 111 (a 2 6 4 — 0462) 
+ (0164 - 0461)111(0263 - a 3 6 2 ) 

Similarly, if we consider the alphabet A = we get 

e (°") a Ml)M2) ' ' ' a V(2n-l)«<T(2Ti)^l ' ' ' ^*2n 

il<---<«2n cre©2n 

= n ( 1 + (°« ~ ^ ) ( is ) 

i \ i<j J 

and the Wick formula gives 

Y e(o-)a CT (i )o .(2) • ■ • a CT (2n-i) CT (2n) = Pf uj (an - a*fc)i<fc,z<2n (19) 

Remark that although (|l~7|) looks like the specialization = a^- of ([19]), 
this is not the case, since this specialization is not a shuffle homomorphism. 
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Example 3.2 with 1 < i,j < 4, the formula reads 

(a i2 - a 2 i) 111 (a 34 - a 43 ) - (ai 3 - a 3 i) 111 (a 2 4 - 042) + (ai4 - 041) 111 (a 23 - a 32 ) 

= Oi 2 a 34 — 012^43 + 034012 — 043 a 12 — 2 i0 3 4 + 2 i04 3 — 34 a 2 i + 04 3 2 i 
— »13a 2 4 + 013^42 — &24013 + 042013 + 3 ia 2 4 — 3 iO4 2 + a 24 a3i — a4 2 a 3 i 

+ ai4«32 — £114^23 + a32«14 — fl23 a 14 — a41«32 + 04i0 23 — 32 a 4 i + 23 04i. 

Similar identities can be obtained by expanding a bosonic version of the 
series (|12]). Let {^} be a set of commuting symbols verifying = 0, and let 
Q be a symmetric matrix with zero diagonal. The series 



7 n(> ^q-ju) 

i \ i<j J 



(20) 



has as coefficients 



UT) = Hf (Q lkH ) (21) 



where I = (ii, - ■ • , im)- Specializing this to the shuffle algebra over variables 
ciij and taking Qm = dki + aik if k ^ I and Qkk = 0, we obtain 

^ a a(l)er(2) ' ' " a CT (2n-l)<j(2n) = Hf m (Q k i) 1 < fc i < 2n ■ (22) 



Example 3.3 Let A = {aij}i<i,j<4- We have 

(a i2 + a 2 i) LU (a 34 + a 43 ) + (a i3 + o 3 i) 111 (a 24 + a 42 ) + (a i4 + 041) LU (023 + 032) 
= ai 2 a 34 + ai 2 a 43 + 034012 + 043012 + a 2 ia 34 + a 2 ia 43 + a 34 a 2 i + 043021 
+ oi 3 a 24 + 013042 + 024013 + a 42 ai 3 + a 3 ia 24 + a 3 ia 42 + o 24 a 3i + a 42 a 3 i 

+ o i4 o 32 + oi 4 a 23 + a 32 ai4 + a 23 ai 4 + a 4 ia 32 + a 4 ia 23 + a 32 a 4 i + 023041- 

3.2 Other Hafnian and Pfaffian identities 

The simplest one in the series is 

^2 ' ' ' a °( 2n ) = (2 n )\\ m m (Q kl ^i<k,i<2n ( 23 ) 

where Qki = ak a i + o-io-k if k ^ I and Qkk — 0. It amounts to the associativity 
and commutativity of the shuffle product. 
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Example 3.4 With A — {a\, a 2 , a 3 , a 4 }, we have 

i [{a\a 2 + a 2 a 1 ) 111 (a 3 a 4 + a 4 a 3 ) + (ai<i 3 + a 3 ai) 111 (a 2 a 4 + a 4 a 2 ) 

+ (a 1 a 4 + a 4 ai) 111 (a 2 a 3 + a 3 a 2 )] = aia 2 a 3 a 4 + aia 2 a4a 3 + aia^a^ + aia?,a4a 2 + 
a\aia 2 ay, + a\a^a 3 a 2 + a 2 a\a^ai + a 2 a\a^aj, + a^zdia^ + a 2 a 3 a 4 o 4 + a 2 a 4 ctia 3 + 
02040301 + a 3 a 2 aia 4 + a 3 a 2 a 4 oi + a 3 a\a 2 ai + a^aia4a 2 + azdACLia 2 + a 3 o 4 a 2 ai + 
o 4 a 2 a 3 a 4 + a 4 a 2 ctia 3 + a 4 a 3 a 2 a 4 + a^a^ai^ + a 4 aia 3 a 2 + a 4 aia 2 a 3 . 

Indeed, the left hand side is just |(3 a\ Lila 2 LLIa 3 LUa 4 ). 

We can find other identities by considering the following generating series 
in Ax[x](^)' wri ere x = {xi, x 2 , • ■ • } is a family of commuting variables: 

T(x, V) = J2 ^ = II t 1 + x ™) ( 24 ) 
/T • 

where / j means that / is an increasing sequence. Since rjf = 0, one can 
write 

T{x,rj) = e H = Y[e x ^ (25) 

i 

where 

H = x iVi + \ ^[xiVi, x jVj] = ^2 Xi7]i + x i x i 1 liVj ( 26 ) 

l j<jr I i<j 

is the Hausdorff series, here terminating in degree 2. On another hand, one 
has 



ti 



= J d m T{x, rf) (27) 



where / = (i n , i n -i, • • • , h) if I = • • • , i n )- If ^(-0 is even (/(/) = 2n), we 
have as well 



ti = J drjj exp I ^ •'•,■' ,//,//, J , 



(2? 



and since 



exp <^ x i x jViVj \ = Yl Yl x ^ x n ■ ■ ■ x in x jnVhVh " " " Vi„Vj n , (29) 



i<j ) n i\<3\,---,in<3n 



8 



J.-G. Luque and J.-Y. Thibon 



we have finally 

ti = Pf {Qki\< kyl < 2n (30) 

where 

Qki = { —Xi k Xij if I < k (31) 
if I = k. 



If l{I) = 2n + l is odd, ([27p gives 



(32) 



But 

e Ei*iW = i + (33) 

i 

and then 

ti = V(-l) fc+1 x lfc Pf (Q y ) i<«< a ^ . (34) 
k 

Let us now specialize this formula to the shuffle algebra. In /\k(a) w (v)^ the 
generating series 

s s ■ ■ ■ a v«)^ = n ( x + Gi ^) ( 3s ) 

gives 

^ JPf us (Qki)i< k ,i< n if n is even, 

2^ a.d)- S E( _ 1)P+ i OpWPf „,(£«)!<*.«<» ifnisodd, W 



where 



afcllla; if k < I 
Qm={ -afcLUa/ if I < k (37) 
if k = I 
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Example 3.5 For A = {a 4 , 02, 03, a 4 }, we have 

(aia 2 + 0201 ) LU (a 3 a 4 + a 4 a 3 ) - (0103 + a 3 ai) LU (a 2 a 4 + a 4 a 2 ) + 
(aia 4 + a 4 ai) LU (a 2 a3 + a 3 a 2 ) = aia 2 a 3 a 4 + aia2a 4 a 3 + aia3a 2 a 4 + aia3a 4 a 2 + aia 4 a 2 a3 + 
aia 4 a3a 2 + a 2 aia3a 4 + a 2 aia 4 a3 + a 2 a3aia 4 + a 2 a3a 4 ai +a 2 a 4 ai03 + a 2 a 4 a3ai + a3<2 2 aia 4 + 
03020401 +03aia 2 a 4 + a3aia 4 a 2 +a3a 4 aia 2 + 03a 4 a 2 ai +a 4 a 2 a3ai + o 4 a 2 ai03 + a 4 a3a 2 ai + 

a 4 030ia 2 + a 4 aia302 + a 4 ai0203 

Clearly, the left hand side reduces to a 4 LUa2 Lila3 LLIa 4 , and the formula does not bring any 
new information. 

For completeness, let us remark that there is a bosonic version of ( p4[) in 
K(A) m [£], where 1T1 denotes the antishuffle product on K(A) (i.e. 1T1 = Lil 9 
for q — — 1, where the g-shume is defined by the recursive formula au\h q bv = 
a{u\n q bv) + q^b(au[n q v)) 

Yl (1 + X&) =J2Yl e ( a ) a V(D ' ' ' a Mn)& (38) 

i it o"GS„ 

leads to the following equality in the antishuffle algebra K m (A) 

V- /n J Hf {Qki)i< k ,i<n if n is even 

^ Z^b a p m Hf hi (Qfci) i<M<n if n is odd 

where = = a fe a/ - a^a fc if fc < Z. 

Expanding the hafnians, we see that this identity amounts to the asso- 
ciativity and anticommutatitity of the antishuffle. 

Example 3.6 Let A={a±, a 2 , 03, a 4 }. We have 

(aia 2 — a 2 ai) in (a3a 4 — a 4 a3) + (0103 — 0301) ITI (a 2 a 4 — a 4 a 2 ) + 
(a 4 a 4 — a 4 a 4 ) ITI (a 2 a3 — a3a 2 ) = aia 2 a3a 4 — aia 2 a 4 a3 — aia3a 2 a 4 + aia3a 4 a 2 + aia 4 a 2 a3 + 
— aia 4 a3a 2 — a 2 ai03a 4 + a 2 0ia 4 03 + a 2 030ia 4 H — a 2 a3a 4 ai — a 2 a 4 aia3 + a 2 a 4 03ai + 
— a3a 2 aia 4 + 03a 2 a 4 ai + 0301020^1 — 030ia 4 a 2 — 030 4 oia 2 + 03a 4 a20i — a 4 a2030i + 

a 4 a 2 aia3 + o 4 a3a 2 ai H a 4 a3<2ia 2 — a 4 aia3a 2 + a 4 aia 2 a3. 

Indeed, the left hand side reduces to 01 ITI o 2 ITI 03 ITI a 4 . 

3.3 Applications to iterated integrals 

Setting ai = bi in ([LTD, we find 

(40) 
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Using Chen's lemma (|J), we recover de Bruijn's formula (|J). Indeed, taking 
as alphabet A = {1, . . . , In) the subscripts of our functions fa, and applying 
the linear form ( ) to both sides of (p0|), we obtain 

e{<T)(v) = (Pi m (ij-ji))=Pf((ij-ji)) (41) 

0"€<32n 

which is exactly (|J). 

Formula ( |TT| ) can also be interpreted in terms of iterated integrals, and 
the corresponding identity seems to be new. Taking 4n functions fa, . . . , fa n , 
■01 j • • • > V^n of 2n variables, we obtain 

det (fa(xj)\ipi(xj)) dx\ ■ ■ ■ dx 2n = Pf (Qij)i<ij<2 n 

(42) 

a<zi<--<a;2ri<6 

where (^(a^)!^^)) denotes the 2n x 2n matrix whose ith row is 
[fa(xi),fa(x 2 ), ... , fa(x 2n -i), ^i(x 2n )], 

and 

Qij = J J lfa(x)ipj(y) ~ fa(x)ipi(y)]dxdy . (43) 

a<x<y<b 




Example 3.7 Let us consider two families of complex numbers (xi)i>i and {yi)i>\ (each 
Xi being non zero). We construct for each i two functions 

f . {t)= i_ t ^-i 



and 



The iterated integral 



= — ■ 

3/i 



/„ = J ' " J det (/ife)l5ife))*l ' ' ' 

0<ti<t 2 <— <t 2 „<l 

can be written as the Pfaffian 

I n = P f 

\ t" «£j 

The calculation of this Pfaffian, due to Sundquist [0~7|], is discussed in section I 
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Application of Chen's lemma to ([19]) gives the second de Bruijn formula 
I), and we can deduce from (|23|) a "de Bruijn like" equality for permanents 

per ((j) i (t j ))dt l ---dt 2 n = 

1 



(2n)!! Hf I J J ^ x ^iv) + <Pj{x)<pi{y)]dxdy | . 

a<x<y<b 

(44) 

which is not very informative, as both sides are clearly equal to the product 
of the integrals f (pi(x)dx. The corresponding identities for determinants 
and permanents of odd order can be obtained in the same way. 

Example 3.8 We can write down a permanental version of Wigner's integral also 
discussed in section || 



/»= /•■■/ peT(e x ^)dt 1 ---dt 2 



a<t 1 <t 2 <-<t 2 n<b 

which evaluates to 

T 1 uf/ft" a i)( b J - a i) \ TT hi ~ <>>i 



(2n)!! \ XiX 
with ai = e aXi and bi = e bxi . 

Example 3.9 The equality 

— \ = (45) 

can be easily shown by induction on m, remarking that 

2 " = Xl + , " + :Em . (46) 

But we can also see it as an integral of type (pi]). Indeed, direct calculation of 
I(xi, ■ ■ ■ ,x m ) = J "J pcr (t{ 3 ^dtx-'-dtm 

0<ti<---<t m <l 

gives 



I(xi, ■ ■ ■ ,x m ) = Y — 



x tr(l) x <j(l) + £<r(2) ^o-(l) + 1" x a(m) 

On another hand, 

/(xi, • • • ,x m ) — Ai ■ • • A m 

with A, = J^t^^dt = i. 
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Finally, applying Chen's lemma to ([22] ) we obtain 
per I V>i dh ■ • ■ dt n 



a<ti<--<t n <b 

= Hf ( I [hix^ix) + faixWiix)] j ( 1.7) 



Example 3.10 The hafnian analog of example 3.7 is 

' Vi + Vj 



I(x, y) = J "J P er (y^T rT) dtl ' ' ' dtn = Hf 

0<tl<— <t„<i 



3>i — f - £Cj 



If we compute I{x\, • • • , X2n, yi, ■ ■ ■ , yin) directly, we find the symmetrization identity 

Va{l)VaCi) ■ •■y<r(2n-l) = Rf f Vi + Vj 

<reS,„ + X <y(2))i x a(l) + X<t(2) + x a(3) + x a(A)) ' ' ' ( x a(l) H 1" £ CT (2r>)) \ X i + X j 

4 Examples 

4.1 Pfaffian identities 

The classical identity of Schur (see |12| for an illuminating discussion of 
this identity and of similar ones) 

•' /.• ~~ ■''/ \ j~j / •' /.• — •' / 

k<l 

(where k, I = 1,... ,2m) has been recently generalized by Sundquist fl7| , 
giving evaluations of 

pu'^zvl) and Pf ( ux-w 



pf 1 = TT ^— ii ( 4 j 



X k + Xi) \l + X k Xi 

and also by Ishikawa and Wakayama JlO | , giving the interpolation 



Pf 



Vk ~ Vl 



a + b(x k + X[) + cx k xi 



Actually, all three identities are equivalent. Indeed, one can obtain Pf 
from Pf f by" setting z k — \ f f^f^ l an d then setting x k — a + /3t 



k 
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one obtains Pf 



yk-yi 



Therefore, to relate these identities 



K (l+a 2 )+aP(t k +ti)+P 2 t k ti 

to the previous considerations, it will be sufficient to consider one of them 
say, the one for Pf 



Vk Vl ). To this aim, we consider the de Bruijn integral 



'2m 



det [e Xktl ] dh 



a<t\<---<t2m<b 

where Q = {Qui) is given in terms of = e aXk and h = e 

1 



dt 2m = Pf(Q) 
by 



(49) 



Q 



kl 



XkXl(x k + Xi) 



Ul Vl 



(50) 



and Uk = b k — cik, Vk = Xk(a,k + The integral ([49"D, which is of interest 
in random matrix theory, can be evaluated directly, without the help of de 
Bruijn's formula. This is a rather tedious calculation, originally due to E. 
Wigner (unpublished), which can be found in Mehta's book |13| (see p. 455, 
A. 23). This calculation results into an expression different from the previous 
one (||13||, Theorem 10.9.1), which in the present notation reads 



'2m 



2m 



detixl'ik + i-iya,)). 



(51) 



Comparing the right hand sides of fl51~|) and ([49]) , and substituting bk — 
bk + a k = x^Vk, one finds 



Ui 




x\u\ 


x\v\ 


xfui 


■ xl m - 2 Vl 


U 2 


v 2 


x\u 2 


x\v 2 


x\u 2 


■ x 2 2 m - 2 v 2 


«2m 


V 2 m 


2 


2 


4 

%2m U 2m 


- h 2m u 2m 



i<j 



UiVj — UjVi 



(52) 



which is one of Sundquist's formulas (when one sets Ui = 1 and Vi = yi). 

In Section |3|, we have interpreted de Bruijn's formula as a combinatorial 
property of shuffles. In the following, we will observe that (at least) parts of 
Wigner's argument can be interpreted similarly. 

To establish ([51]), Wigner starts with the simplest integral 



Jn {X\ , . . . , X n ^ 



Mxi-{ HnXr 



dt\ • • • dt n , 



(53) 



a<tl<-<t2m<b 
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which yields I n by antisymmetrization: 

I n (xi, ■ ■ ■ ,X n ) = AJ n (xi, . . . ,X n ) = s(a)J n (x a (X),. ■ ■ , Xg(n)) ■ (54) 



If one computes the first integrals J\, J 2 , ■ ■ ■ , an induction pattern emerges, 
and one arrives at the expression 

Jn\XXi ■ ■ ■ ; Xji) 
n 

^2(-l) k a k ■ ■ ■ a 2 ai • b k+1 ■ ■ ■ b n R(x k , x k -i, ■ ■ ■ , x 1 )R(x k +i, x k+2 , ... ,x n ) 



k=0 



where we have set 



(55) 



R(z u . . . ,z r ) = — ■ r j — ■ ■ -. (56) 

z\\z\ + z 2 ) ■ ■ ■ {zi + z 2 H z r ) 

The rest of the calculation relies upon the two identities 

n 

^2(-l) h R(x k , x k -i, ■■■ , xi)R(xk+i, x k+2 , ... ,x n ) = (57) 

k=0 

and 



n 1 — 

AR(x 1 ,...,x n )=T[-T[^-^. (5f 



1=1 1<J J 



Remark that if we set <j> z {t) = t z , we can write 



R(z 1: . . . ,z r ) = (zi ■ ■ -Zr) = J '"J 0*i (*i) • • • 4zr(tr)dti ■ ■ • dt r (59) 

0<*i<-<t r <l 

so that for u = x k ■ ■ ■ x±, v = x k+ i ■ ■ ■ x n (regarded as words), R(u)R(v) = 
(u)(v) = (uWv). 



With this at hand, it is clear that fl57|) amounts to the fact that the map 



S(w) = (— 1)^vj, where w is the mirror image of the word w, is the antipode 
of the shuffle Hopf algebra. This means that in the convolution algebra of 
Km (A), i.e., End gr (i^ m (A)) endowed with the product f*g = sho {f®g) oC, 
where sh(u <g> v) — uuiv and C{w) = J2 U v=w u ® v i one nas 

S*I = I*S = Loe 
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where i is the unit and e the counit, so that S * I of a polynomial is equal to 
its constant term. Hence, for any non empty word w, 

(-l) H wWv = 0, 

uv=w 

which yields Equation ( |5*gD can be obtained by writing 

AR(x h ... ,x n )= J ■■■ J det^fo))^ . . . dt n = Pf (A) 

0<ii<---<t„<l 

where 



m 3 = [ dtj f' <»,ir; 1 V ') 

Jo </o 
and applying Schur's identity 



There is also an interesting interpretation of (|52|), involving the hype- 
roctahedral group. Let us introduce a family of functions labelled by the 
alphabet A = {n, ... , 2, 1, 1, 2, . . . , n} 

4> k {t) = b^ 1 and (j)- k (t) = -a^ 1 , (60) 

so that (|55|) can be written as 

J n (xi,... ,x n ) = (j^k-- ■21\llk+ l---n\ (61) 

\fc=0 / 

and 

A7„ = / e(a)o-^fc---2lLUfc+l---n\ (62) 
\o-Ge„ fc=0 / 

If the word w occurs in k • • • 21 111 k + 1 • • • n with the letters k, . . . , 1 
at positions it,.. . ,ik, the sign of the underlying permutation a, obtained 
by erasing the bars, is e(a) = ( — 1) " z " , where ||/|| = — 1)- We can 

furthermore interpret such a word w as an element of the hyperoctahedral 
group B n , generated by permutations and by the involutions Tj : i \— > i and 
T i{j) = 3 f° r 3 7^ i- The double sum in (|62]) can now be interpreted as an 
element of the group algebra ZS n , and since Aa = e(a)A, we have for w as 
above 

Aw = (-l)^Ar h ---T ik . (63) 
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If we extend the action of B n to polynomials in x^y^a^ bi by setting Xi = X{ 
Vi — Hi and bi = — a*, (|62D can now be rewritten 



AJ n = Ao(l+T 1 )o(l-r 2 )o...o(l-r 2 " Wh-'-hn, 



Xi(xi + x 2 ) ■ ■ ■ (x 1 H h x 2m ) 

(64) 



that is, 



where 



AJ n = A(QR) (65) 



2m 



Q = H(b i + (-l) i a i ) (66) 



i=i 



Let x p = a??^ • ■ ■ x™ 1 , so that A = *4.(.x p ). Then, we see that fl52| ) is 
equivalent to the symmetrization identity 

A{RQ)A{x p ) = A{R)A{Qx p ) . (67) 



4.2 Anticommutative symmetric functions and altern- 
ing quasi- symmetric functions 

The fermionic Wick formula admits an interpretation in the language of non- 
commutative symmetric functions ||. Let us recall that the starting point 
of this theory is the observation that the classical algebra of symmetric func- 
tions is a polynomial algebra Sym = C[h±, h 2 ■ ■ ■ ] in an infinite number of 
indeterminates (here, the hi are the complete homogeneous symmetric func- 
tions). One defines the algebra Sym of noncommutative symmetric functions 
as the free associative algebra C(Si, S 2 , • ■ ■ ) over an infinite sequence of let- 
ters, together with a "commutative image" homomorphism S n — ► h n . In 
other words, Sym is the tensor algebra of the vector space V = CS n and 
Sym is its symmetric algebra. 

Here, we propose to have a look at its exterior algebra, which we will 
call ASym (anticommutative symmetric functions). It is the image of Sym 
under the specialization S n — > S n (A) = r\ n (that is, the letter A is used here 
as a label for the image F(A) of any noncommutative symmetric function F 
by the ring homomorphism S n i— > r) n ). 

The dimension of the homogeneous component of degree n of ASym is 
equal to the number of partitions of n into distinct parts. As an algebra, 
ASym is rather trivial, but the dual Hopf algebra looks more interesting. To 
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study it, we can follow the same strategy as in |J, and start with the Cauchy 
kernel. In the notation of |H| 



a,(XA) = ]Ja Xk (A) = U £4X 

k>l k>li k >0 

= E E ^eWM^XjU,-^ (68) 

r ii< — <i r \cr£&r / 

= E E (e^) m ^ x )V( a ) 

where the Mi(X) are the quasimonomial quasi-symmetric functions, and la 
means • • • , V(r))- The coefficients 

Vj(X) =J2<*WUX) = I dma^XA) (69) 

a J 

can be evaluated by means of the fermionic Wick formula. Introducing the 
operators 

d f 

4>k = ^— = drj k and <f>* k (f) = rj k f 

07] k J 

one has 

V I =(Q\(j> ir ---<f>iMXA)) (70) 
where |0) = 1 and a x {XA) = JJe* (a!i) |0) if a x (A) = (here the generat- 

i>l 

ing series $(x) of the second kind power sums reduces to xr\\ + x 2 n 2 + •••). 
The Hausdorff series for the product of exponentials is now 

F = $(Z 2 ),---) 

fc>l i<j 



£ M n (x)0; + ~J2 ( M ^( x ) - 

5 + (7i) 



2 

n>l fc^i 



2 

fc^Z n>l 



where Q M = M fci - M ifc . Hence, 

Vi(X) = J d m e H = (0|e E «>° M ™^e^*W*| ) (72) 
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since the quadratic terms are central, so, expanding the exponential e^ n Mn ^n^ 
we obtain the following expression for the "alternating quasisymmetric func- 
tions" . 

Vr(X) = <a)M Ia (X) = (<j> ir • • • <j> h ) 

(7 

Pf (Q ikil )i<k,i<2m if r = 2m is even 

.^(-^(^"A"^) ifr = 2m + lisodd. 

This formula is similar to Schur's expression of his symmetric functions 
Qx as Pfaffians of the Quj)- Remark that these symmetric functions are also 
labelled by strict partitions, and that Schur's Pfaffian can also be derived 



from the fermionic Wick formula 19 . 



5 Generalizations 
5.1 Hyperpfaffians 

The Pfaffian of a skew symmetric matrix of even order can be expanded as 

Pf ( M ^l<i,j<2n = E <°) M <r(lW) ■ ■ ■ ^(2n-l M 2n) (74) 

where the sum is over 



e 2 n,2 = {ae e 2n \a(2i + 1) < o(2i + 2), a(2p - 1) < a{2p + 1), 

0<i<n-l,l<p<n-l} (75) 

We can extend this definition by setting 

<£fcn,fc = {a e & kn \a(ki + 1) < ■ • • < a(k(i + 1)), 

a(k(p - 1) + 1) < a(kp + 1), < i < n - 1, 1 < p < n - 1}. 

If we consider variables M h ... ik verifying M ia{1) ... itj{k) = e(a)M h ... ik , the poly- 
nomial 

Pf [fc] (M n ... J J 1 < nj ^ 4fc < nfc = <^)Ma(l).^)--- M -((n-l)k + l)...,(nk) (76) 

appears as a natural generalization of the Pfaffian. This polynomial has been 
called the hyperpfaffian of the alternating tensor M by Barvinok JTJ . 
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Similarly, higher order Hafnians, or hyperhafnians, can be defined by 

Hf ^(M il ...jji< ilj ... jifc < rlfc = ^2 Mr(l)...<r(fc) • • •Mr((n-l)fc+l)...<7(nfc)- (77) 

Alternatively, these definitions can be presented as follows. Given an 
alternating tensor M of order k as above, consider the element 

Q M = Yl M h-i k Vh ■ ■ ■ Vh (78) 

l<ii<...<ik<nk 

of the Grassmann algebra of rank nk. Then, in this algebra 

n^ = n\PfW(M)r ]lV2 --- Vkn . (79) 

When k = 2, this is one of the classical definitions of the Pfaflian. 

Similarly, when M is a symmetric tensor of order k, again over an nk 
dimensional space, consider the element 

G M = J2 M n-iu^--<i* ( 8 °) 

l<«i<...<ifc<nfc 

wher £i, . . . , £ fcn are commuting symbols satisfying = 0. Then 

G n M = nm&\M)^ 2 ---i kn . (81) 

5.2 Composition of hyperpfaffians 

The aim of this subsection is to compute the hyperfaffian 

Pf [ ] (Pf (Oj ,j)i,je{k u - Mm}) kl ... ,fe 2m6 {l,... ,2mn} 

where A = (aij)i<ij<2 mn is a skew symmetric matrix and m, n > 0. In order 
to express this as a classical Pfaflian, we consider the 2— form 

U = a ijViVj ■ 

l<i,j<2mn 

The mnth-power of u is by definition of Pfaffians 

w™ = (mra)!Pf (a i j)i<i ij < 2m „ . 

But, on another hand 

uj mn = (m\0) n 
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where 

m 

Remarking that Q n is given by a hyperpfaffian 

Q n = n ,p f [2m] (p f ( Qi .) . . 6{feii ... fc2m} ) ^ ^ . . . ^ , 

by identification, we find 

pf [2m] (Pf (a iJ )j )je { fcl ,... J fe 2m }) 1 < fei ... jjfc2ro < 2TOn = / n A Pf (Oij)l<ij<2mn 



(82) 



Example 5.1 If we set m = 2 and a^j = , we can generalize Schur's formula to the 
4-hyperpfaffian 

pf[4] f n x ^t) =(2n)!! n ^ (*» 

V Ks t<4 ks kt I Ki i<4n 1 3 

5.3 Hyperpfaffian of a sum 



In |[L6|| ) J. Stembridge gives a formula allowing to compute the Pfaffian of a 
sum of 2s x 2s skew symmetric matrices 

s 

Pf{A + B) = J2 (- 1 )' 7h * Pf (Aj)Pf(.Bj<0 (84) 

t=o /=(*!,... ,i 2t )T 

where 1 1\ — i\ + • • • and J c = {1 • • • 2m} — J. 

This identity appears to be a special case of a more general one for hyperp- 
faffians 

n 

Pf [2m] {A + B) = (-l) |/hfcm Pf [2m] (A/)Pf [2m] (S 7 c) (85) 

k =0 I={h<-<i2km} 

where A = (cLi lr - ,i 2n ) and B = (&i lr .- ,j 2n ) are two skew symmetric tensors of 
order 2m over a space of dimension 2mn. 
We set 

^= ( a n-W + & n-w)%---% m - 

1<«1 <--<i 2m <2mn 

The nth-power of a; is 

a;« = n!Pf ^(A + B)t& • • • % mn . 
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On the other hand, we have 

n 
fc=0 

where u A = £ a iv .. i2m r] h ■ ■ ■ n i2m and u; B = £ b h-i 2m Vh • • • Vi 2m - 
Since 

il<-"<t2fcm 

we find 

u," = n! E E Pf [2ml ( A/ ) Pf ^ (Brimi' > 

fc J={*l<-<*2fcm} 

where J c = {1, • • • , 2mn} — J. Finally, 

w ™ = n!^ E (_l)m-*mp f [2m]^p f PM^)^ 

IC{l"-2mn} 

which by identification leads to (|85|) . 

5.4 Minor summation formula 

In [H], M. Ishikawa and M. Wakayama prove a minor summation formula for 
Pfaffians 

E pf ( B tt) det (rfcXJ = Pf (Q) , (86) 

l<fci<---<fc m <n 

where A 1 / '"t™ denotes the submatrix of ^4 which consists of the rows ii • - • i m 
and the columns ji, - ■ • , j' m , 5 is a skew symmetric matrix, T is an arbitrary 
matrix and Q = TB l T. We shall now extend this formula to hyperpfaffians. 

Let A = (ai 1 ...i 2m )i<i 1 ,...,i 2m <2mn be a skew symmetric tensor of order 2m 
over a space of dimension 2mn and T = (Tij)x<i,j<2mn be a square matrix. We 
want to find a skew symmetric tensor Q = (Qi 1 ,...,i 2m )i<i 1 ...i 2 „<2mn verifying 

E pf [2m] (^-w) **(T£XJ = pf [2m] W) > ( 87 ) 

l<fci<-<J:2 m i<2mn 

where -A^...^,^ = (aii,-,i 2TO )ii,-,i 2ro 6{fei,-,fe 2mt }- 

In order to compute such a tensor, we need first to remark that 

a* e«i) = e *o det • • ■ ^ {<^t:i 2mn) ) 

(88) 
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which can be obtained by considering uj k = Ylii=i TikVi an d computing uj kl ■ ■ ■ u kt . 
Now, let us consider 

l<fci<'-<fc 2m <2mn 

and the polynomial 

E ( E «fe--^^(^^£:) 

l<U<-- <i2m<2mt \fci< --<fc 2m / 



A simple calculation gives 
On another hand 



m(t — 1)+1 '"*2mi Vil ^]i2nt 

<76<£2m,2mt 

^ ; ^ ; a ka(l)---k lT (2m) ' ' ' a ka-(2m(t-l) + l)--- a <T(2mt) 

Kfcl<-<fc2mi<2mti CTlGl£2m,2mt 

V e(<7 2 ) det f T r (1) -? (2m) V • ■ det f i«<Mt-i)+i>-^<**r 
/ j \ *j ^ St 1 )" ■'li 2 ™) / V S( 2m ('-i)+i)"'S( 2ml ) , 

X 7/1 • • • r/ 2m4 . (89) 



(72 6 £2 



By formula fl38f), we find 



- L ' ~"~ afe CT (l)--- fe CT (2m) ' ' ' afc «T(2m(t-l) + l)---«<7(2r ■•■/ , 



xdet^--^)^.-.^ (90) 



so that 



= t! E Pf [2ml CV~ det ( T tXJVl ■ ■ ■ V2mt , 

l<fcl<'--<fc2mt<2mn 



which leads to the identity fl57p. 
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5.5 Generalized de Bruijn formulas 

We consider now a commutative ring K and for 2k a non zero even integer, 
a set of scalars {Qi 1 --i 2k }i 1 < -<i 2k - We shall first calculate the coefficients Pj 
of the /\ K (^-series 

f = J2pivi = U[ 1+ E Qh-i 2k Vh---Vi 2k )- (91) 

J| ii \ ii<-<j fc / 

As ?7j 2 = 0, one has 

(92) 

n i 2fc(p-i)+i<-<'2fcp. 1 ^P^» 

i 2fe(q-l) + l <i 2fcq + l. 1 <9<™-l 

and if I = {ii, ■ ■ ■ , %2nk\ is an increasing sequence, the coefficient Pi of r\i is 

ft = Pf M (fl^) • (93) 
Similarly, in the commutative algebra K[£\ ( £f = 0) 



e Hf[2fei K-^)^=n( i+ e ■■<,.) 

J T h V »i<— <»2/t / 

l(I)=2nk 



(94) 



In the remainder of this section, the ring K is the shuffle algebra R{A) m 
for the alphabet A = {aj 1 ...j 2 Ji 1 ,...,i tt >i.We consider here the series 

T{A\rf) = n( 1+ E ( E e ( c7 ) a V(i)-V(2 fe) ) Vh---Vi 2k 

h V ii<— <ifc \cre6 2 fe / 

= E p ^ 

IT 

i(/)=2fcn 

If we set S(ii, • • • , i 2k ) = E CTG e 2fc e(cr)a i£T(1) ... i(T(2fc) , we obtain for each increas- 
ing sequence I = (ii, • • • , i 2fen ) 

Pi — J d m T(A\ V ) 

n<---<hk 

= ^2 e (°")^(^(l)' ' ' ' >M2fc))LU ' ' ' LUS(V(2(„_l)£; + l), • • • ,V(2nfc)) 

0"e<£ 2 n,2nfc 

~~ ^ ] e ( Cr ) a i<T(l)---* CT (2fc) ' ' ' a M2(n-l)fc + l)-"V(2nfc) 1 
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On the other hand, we find by fl93|) 

e^K.tD-M^) I ( 95 ) 



By Chen's lemma, this equality leads a generalized de Bruijn formula. We 
consider a set of 4/c 2 n functions $ = {4>ij} i<;<2fc . Equation fl95|) gives 

1<J <2kn 

2U$(a,6) = / ••• / det(0ij(a;j)| • • • |02fci(^-)) i<i<2nfc rfxi . . . (ix n 

a<a;i<---<x n <6 

V 7 " 1 <xes 2fe 



a<a;i<-<x n <6 



&2k ' l<u,---,i 2fc <2nfc 

(96) 



Similarly there is a permanent version of this identity 
971$ (a, 6) = /■•■/ per(0i i (x : ,)| ■ • • |0 2 fci(^i))i<i<2nfc(ia;i . . .dx n 



a<xi <---<x n <b 



Hf[2fc] / Y, <i>u„ ( u( x )--- ( i ) 2^ (2k) (x)dx 



W(r(2fc) 1 



1 <ii ■ ■ ,i 2 fc<2nfc 

(97) 



5.6 Sums of even powers of the Vandermonde 

We can use formula fl9~ED to obtain the average of certain determinants. For 
example, let us consider n real numbers xi, ■ • ■ ,x n chosen at random from iV 
values yi, ■ ■ ■ with uniform probability. The average of the polynomial 
) (for m G N) can be expressed as an integral 

A(n,m) = J ■■■ J A 2m ( Xl ,--- ,x n )f[ (^25( Xi - y p )j d Xi 

where 8 denotes the Dirac distribution. We can transform A(n, m) into the 
interated integral of single determinant 

A ^ m ) = jfi J - j det (fl(xM{*i)\ ■ ■ ■ \f? m (xj)) dsn ■ --dx n 

— oo<xi<---<x n <oo 
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where 

i-1 



x i-(s-l)n-l tf(s-i) n + i<i< sn 
otherwise 



fi (x \ = { x Ep=i K x ~Vp) if i < " 
Ji{ ' \ otherwise 

and 

/• (*) = 

if 1 < s < 2m. Using formula (|9^) we find 

An,m) = ^f [2m] (Mi 1 .., 2 J (98) 

where if i\ < ■ ■ ■ < i 2 m 

r°° s(h - io \s^n r / \ 7 if for each s < 2m, 
, x s[ 11 ' > i2mj V„, 8(x - y v )dx , 1N . . 

3 otherwise 

(99) 

with 

2m 

• • • , km) = ~ m(2n(m - l) + 2) 

k=i 

and 

The nonzero entries of the hyperpfaffian are therefore the power-sums 

N 

p=i 

with »!<•••< and (s — l)n < i s < sn for each s < 2m. 
In the case where m = 1, we find 

^ 1 ) = ^det^^. (101) 

More generally, if we consider the measure 

*-6&)n(£^)*. 
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we can generalize formula fl98|) to evaluate 



— pf [2m] r M (k) 



where 



At 



(k) 



AT s(n,---,l2m) 

E E 

p=l t=o 



yp (2ky 



t 



(102) 



(103) 



where H t denotes a Hermite polynomial. We can verify that lim A(k, n, m) = 

fc^oo 

A(n, m), as expected. 

Finally, one can generalize formula (BSj) by considering a matrix of func- 
tions 



<2nm 
KK2m 



and the determinant 

P J p(s 1 ,...,s B ) = det {fl(x j )\---\fl m (x j )). 

Its average, when the Xi run over the set • • • , 2/jv} endowed with the 
uniform probability, can be computed as above: 



A(n,m) = — — 



Pf(xi 



) " " ' ; X n 



-00<X\<---<X2m<00 



n / N \ 

II ( ^2 s ( Xi -yp) ) 

i=l \p=l / 



n! 



Pf'(x\, • • ■ , x n )dxi ■ ■ ■ dx n 



-00<Xl<---<X2m<00 



(104) 



where 



\ Sp=l yp)f2mn( X ) fimn( X ) 

By formula we find 



f? m (x) \ 

fimn ( x ) J 



(105) 



A(n,m) = ^P^(M tl ..., i2m ) 



(106) 
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where 

J -°°ae&2m P=l 

N 

To conclude, let us remark that the hyperpfaffians (^Sj) and ( |102| ) can 
alternatively be expressed as hyperdeterminants in the sense of |l[ . We note 
also that in dimension nk, the hyperdeterminant of an alternating tensor A 
of (even) rank k satisfies DET(A) = (Pf *- k \A)) k , generalizing the classical 
case k = 2. Such relations will be investigated in a forthcoming paper. 
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